REDUCED BASIS A POSTERIORI ERROR BOUNDS FOR THE 
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Abstract. We present reduced basis approximations and rigorous a posteriori error bounds for 
the instationary Stokes equations. We shall discuss both a method based on the standard formulation 
as well as a method based on a penalty approach, which combine approaches developed in 0[9] and 
[6] with current reduced basis techniques for parabolic problems 11. 12 15] . The analysis then shows 
how time integration affects the development of reduced basis a posteriori error bounds as well as the 
construction of computationally efficient reduced basis approximation spaces. To demonstrate their 
j^JQf performance in practice, the methods are applied to a Stokes flow in a two-dimensional microchannel 

with a parametrized rectangular obstacle; evolution in time is induced by a time-dependent velocity 
profile on the inflow boundary. Numerical results illustrate (i) the rapid convergence of reduced 
basis approximations, (ii) the performance of a posteriori error bounds with respect to sharpness, 
and (iii) computational efficiency. 
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Introduction. Designed for the real-time and many-query context of parameter 
estimation, optimization, and control, the reduced basis method permits the efficient 
yet reliable approximation of input-output relationships induced by parametrized par- 
tial differential equations. The essential ingredients are: (i) dimension reduction, 
through Galerkin projection onto a low-dimensional reduced basis space; (ii) cer- 
tainty, through rigorous a posteriori bounds for the errors in the reduced basis ap- 
proximations; (iii) computational efficiency, through an Offline-Online computational 
strategy; and (iv) effectiveness, through a greedy sampling approach. 

In this work, we provide an outlook on how reduced basis techniques presented 
in [SI [5] for parametrized saddle point problems may be extended to the time- 
dependent setting. To this end, we consider the instationary Stokes equations. We 
shall discuss both a method based on the standard formulation as well as a method 
based on a penalty approach (see also [8] for initial results) , which combine approaches 
developed in [Jj [9] and [6] with current reduced basis techniques for parabolic problems 
(see, e.g., [ITJ [TH [T5] and also [TSJ[TS]). The analysis then shows how time integration 
affects the development of reduced basis a posteriori error bounds as well as the 
construction of computationally efficient reduced basis approximation spaces. 

Starting from the standard mixed formulation of the instationary Stokes equa- 
tions, we develop rigorous a posteriori error bounds for the reduced basis velocity ap- 
proximations. As in the stationary case presented in [7JI5], they involve the (Online-) 
estimation of coercivity, continuity, and inf-sup stability constants associated with 
the diffusion term and incompressibility constraint; in addition, they now also de- 
pend on continuity constants associated with the mass term. Employing a penalty 
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formulation, we obtain rigorous upper bounds for the errors in both the velocity and 
pressure approximations. As in the stationary case presented in [BJ, they are compu- 
tationally very efficient since they do not involve the estimation of inf-sup constants 
but only depend on coercivity constants associated with the diffusion and penalty 
terms; however, they again also depend on the penalty parameter such that associ- 
ated effectivities increase as we approach the non-penalized problem. To construct 
efficient reduced basis approximation spaces, we consider a POD greedy procedure 
(see [TTJ [TJJ [TS]) that is coupled with adaptive stabilization techniques developed 
in [7]. To demonstrate their performance in practice, the methods are then applied 
to a Stokes flow in a parametrized domain where evolution in time is induced by a 
time-dependent velocity profile on the inflow boundary. 

This paper is organized as follows: In fjl] we introduce the general problem for- 
mulation and its "truth" approximation. We start from a time-discrete framework 
already that allows us to directly recover the settings discussed in [7J [9] and [6J ; now, 
we have a saddle point problem associated with each time step. The time discretiza- 
tion scheme is here given by a backward Euler method. Section [5] then describes our 
reduced basis method. In jj2.1i we define our reduced basis approximation as the 
Galerkin projection onto a low-dimensional reduced basis approximation space. We 
develop rigorous a posteriori error bounds in H2.2I Both reduced basis approxima- 
tions and error bounds can be computed Online-efficiently as summarized in £)2.3I 
This enables us to engage adaptive sampling processes for constructing computation- 
ally efficient reduced basis approximation spaces, which shall be outlined in £12.41 In 
Sj3j we introduce our instationary Stokes model problem. Numerical results in Sj4]then 
illustrate (i) the rapid convergence of reduced basis approximations, (ii) the perfor- 
mance of a posteriori error bounds with respect to sharpness, and (iii) computational 
efficiency. Finally, in fj5j we give some concluding remarks. 



1. General Problem Statement. 



1.1. Formulation. Let X c and Y c be two Hilbert spaces with inner products 
(•lOxja and associated norms || • || Xe = y/(-,-)x c , II ■ Ik = VTvhT: respec- 

tively!^ We define the product space Z c = X c x Y c , with inner product (-,-)z c = 
(•, -)x + (•, -)y b and norm || • \\z c = \J (•, -)z B - The associated dual spaces are denoted 
by X' e , Y:, and Z' e . 

Furthermore, let T> C K™ be a prescribed n-dimensional, compact parameter set. 
For any parameter fj, € T>, we then consider the continuous bilinear forms m(-, •; /i) : 
X c x X c ^ R, a(-, •; /x) : X c x X c ->• M, and b{-, •; n) : X c x Y c ^ R0 

n , x m(u, v, u) _ . . 

^(p) = sup sup ,l ~> < oo, V ^ € V, 1.1 

uex c vex c \\u\\x a \\v\\x e 

7 «(/i) = sup sup ^TT < °°' V » e P ' ( L2 ) 
u<£X a vex \\u\\x e \\v\\x e 

= sup sup { b ^ff < oo, V m € V, (1.3) 



1 Here and in the following, the subscript e denotes "exact". 

2 For clarity of exposition, we suppress the obvious requirement of nonzero elements in the de- 
nominators. 
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as well as c(-, ■; /i) : Y c x Y c —> R, 

7c e ( M ) = sup sup iT^TTT <w > V ^ e2? - (1-4) 
pen gen INInlklln 

We moreover assume that a(-, •; /i) and c(-, •; /i) are coercive on X c and Y c , respectively, 
aSMS «S.l^ >0 ' (L5) 

a ^^^W >0 ' VMe2? ' (L6) 

m(-, -;^) is symmetric and positive definite, 

m(v,v;n) > 0, V0^el c , V/i€P, (1.7) 
and •; /i) satisfies the inf-sup condition 

/3 C ( M ) = inf sup n-^^TT > 0, V n G V. (1.8) 

By (|1.2[) . (jl.5[) and (|1.4p . (|1.6[) . the bilinear forms a(-, -;/x) and c(-, t nen provide 
with || • \\x e ,n = and II • \W B ,n = y/c(-,-;fi) energy norms on X c and F e , 

respectively, which are equivalent to || • ||x e and || • ||y e for any fi G T>. Furthermore, as a 
symmetric and positive definite bilinear form, m(-, •; fj,) defines an inner product on X c 
for any parameter /i £ T>; the associated norm shall be denoted by || • || M = y/ m(-, •; fi). 

We further assume that we are given a time interval [0,T], T > 0, and linear 
functionals /(•; /j,) G C°(0, T; and #(•; fi) G C°(0, T; F e ') for all /i G D ; for a vector 
space V, C°(0, T; V) here denotes the space of V- valued functions of class C° with 
respect to t G [0,T]. Throughout this work, we directly consider a time-discrete 
framework: We divide the time interval [0, T] into K subintervals of equal length 
At = T/K, and define t k = kAt for all k = 0, . . . , K ; for notational convenience, we 
also introduce K = {1,...,K} and K = KU{0}. We then set / fc (-; fi) = f{t k ;n) G X£ 
and /i) = ,g(i fe ; fi) G F e ' for all & G K , /i G V. 

For e > 0, we now consider the following "exact" — more precisely, semi-discrete 
- problem resulting from a backward Euler method (see, e.g., [H Q3J [23l ES]): For 
any given parameter fi & V,we find u £ - k (/i) G X c and p £ - k (fi) £ 7 C , i € K, such that 
u|'°(/i) = C| and 

+ a(u £ /(p),v;n)+b(v, P £ c k (ti);n) = f k {v;fi), V v £ X c , 

keK. (1.9) 

b{u £ e ' k (ti),q;ti) -ec(p e c < k (n),q;n) = g k (q]n), V q G F e , 

Even though we here use a common notation for simplicity in exposition, we point 
out that (|1.9p states very different problems for s = and e > 0, respectively. For 
e = 0, we also denote = u°' fc (^i), k G K , and = Pe' fc (/-0) k £ K, 

for all /i G 2?. For e > 0, corresponding to our discussions in [6], (| 1 .9[) can be 
considered as a perturbed or regularized version of the problem associated with e = 0; 



3 We here assume zero initial conditions for simplicity; note that nonzero initial conditions can 
be handled as well without much difficulty (see I12| ). 
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in this case, we therefore call (ul ,k (/i) , pl ,k (fi)) , k e K, also the penalty solution. 
Since these problems differ considerably in their general nature (cf. [7J and 0), we 
shall often treat them separately in the following analysis and explicitly distinguish 
between the two cases e = and e > 0. From (|1.5[) and (|1.7[) . the bilinear form 
jjm(-, •;//) + a(-, •;//) is coercive on X c for any fi £ T>. The problem (| 1 .9[) is thus 
uniquely solvable for (u k (fi),p k (fj,)), k € K, and (ul' k ([i) , pl' k (/*)), fc € K, as a saddle 
point problem according to [7] and [B] , respectively. 

1.2. Truth Approximation. We now introduce a high-fidelity "truth" approx- 
imation upon which our reduced basis approximation will subsequently be built. To 
this end, let X and Y denote finite-dimensional subspaces of X c and Y c , respectively. 
We define the product space Z = X x Y , and denote by AT the dimension of Z. We 
emphasize that the dimension Af is typically very large. These "truth" approximation 
subspaces inherit the inner products and norms of the exact spaces: (•, -)x = (■, Ox , 
|| ■ \\ x = || ■ II*., (-, -)y = (•, On, II • Ik = II • Ike, and (•, -)z = (•, -)z., \\ ■ \\z = \\ ■ \\zl 

Clearly, the continuity properties (11.11) . (|1.2[) . (|1.3p . and (|1.4p are passed on to 
the "truth" approximation spaces, 

7m (m) = sup sup ' ' < oo, V/i€D, (1-10) 
uex^ex ||u||x|M|x 

. . a{u, v;u) . . t~\ i i\ 

7 a (/i) = sup sup •jj— J, 7—r— < oo, V /i G T>, (l-H) 

uex uex |p||x|M|x 

7b(^j = sup sup t— r — j— j— < oo, VfieV, (1-12) 
gey »ex IMMMIx 

7c (/i) = sup sup i^TlT <(X3 ' V M e2?5 (1.13) 
pexqex \\p\\Y\\q\\Y 

so are the coercivity properties (| 1 . 5|) and (|1.6p . 



a(tvu^i) 

eeX ||»|& 



a a (/x)= mf V |,;,||2 >0. V/ieP, (1.14) 



a c ( M ) = inf C ^'^ M) > 0, V/ieP, (1.15) 
gex ||g|| y 

as well as the inner product m(-, •; /x) and associated norm || ■ || M , 

m(v,v;iJ,)>0, VO^uel, V/ieD. (1.16) 

Thus, || • Hx/i = II ' \\x e ,fj, and || • ||y. M . = || • ||y clAt define norms on X and Y, respectively, 
which are equivalent to \\-\\x and || • ||y for any fi ET>. Wc now further assume that the 
approximation spaces X and Y are chosen such that they satisfy the Ladyzhenskaya- 
Babuska-Brezzi (LBB) inf-sup condition (see, e.g., [2]) 

P(ji)= inf sup >0°Qi)>O, V/.6D, (1.17) 

gey^ex ||?||y|M|x 

where /3 (/i) is a constant independent of the dimension J\f. 

Our high-fidelity "truth" discretization for (|1.9j) now reads as follows: For e > 
and any given // € P, we find u e ' k (fi) € X and p e ' k (n) £ Y, k £ K, such that 
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u e '°(/i) = and 

+ a(u e < k (ti),v;n) + b{v, P £ ' k (ii);ti)= f k (v ;t i), Vwel, 

, , , keK. (1.18) 

b(u e ' k (iJ,),q;iJ,) -ec(p E ' k (n),q;n) =g k (q;n), V q G Y, 

In case of e = 0, we also denote u k (p) = u°' k (fi), k £ K , and p fe (^) = p°' k (/J.), 
k £ K. As the exact problem in £11.11 the problem (|1.18[) is uniquely solvable for 
(u k (p) , p k (//)) , fc e K, and (u^ fc (^),p e ^(^)), fc £ K, according to [7] and [6], respec- 
tively. 

2. The Reduced Basis Method. We now turn to the reduced basis (RB) 
method, discussing the approximation procedure, rigorous a posteriori error estima- 
tors, and the construction of stable approximation spaces that capture the causality 
associated with the parameter dependence as well as with evolution in time. 

2.1. Galerkin Projection. Suppose that we are given a set of nested, low- 
dimensional reduced basis approximation subspaces X^ C X^+i C X and Yv C 
Y/v+i C Y, N £ N max = {1, . . . , iV max }. We denote by N x and N Y the dimensions 
of Xn and Yv, respectively, and the total dimension of Zm = Xn x Y/v by A^- = 
Nx + Ny- The subspaces Xn, Yv, and Z/v again inherit all inner products and 
norms of X, Y, and Z, respectively. The reduced basis approximation is then defined 
as the Galerkin projection with respect to the truth problem (|1.18|) onto these low- 
dimensional subspaces: For e > and any given jj, £ Z>, we find u 6 ^ k (p,) £ Xn and 
PiV 0-0 G Yv, k £ K, such that u^°(/i) — and 

+ a(u% k (p,),v N ;fi) + b(v N ,PN k (n);n) 
b(u% k (fi) , q N ; - £ c(p^ fc (p),q N ;p.) 

Again, we denote u k N (p) = u^{p), k £ Ko, and p%(p.) = (fi), k £ K. 

The discrete reduced basis system now essentially behaves as in the stationary 
case: We recall (see [7]) that a pair (Xjv, Yv) of reduced basis approximation spaces 
is called stable if it satisfies the inf-sup condition 

P N (p) = inf sup ^Tif^ >0, V M eD. (2.2) 
wGYjv VN ^x N HQ'jvIIxII'WjvIIx 

In case of e = 0, (|2.ip is then uniquely solvable for (u%(fi),p%(fi)), k £ K, if and only 
if the reduced basis approximation spaces X^,Y^ are stable; in case of e > 0, (|2.1[) 
is uniquely solvable for (u e ^ k (p) , p e ^ k (fi)) , k £ K, for any choice of Xjv, Yv (see [2 [6]). 

2.2. A Posteriori Error Estimation. We now develop upper bounds for the 
errors in our reduced basis approximations that are rigorous, sharp, and compu- 
tationally efficient. The errors shall be measured relative to the respective truth 
approximations. 

In this section, we assume that the low-dimensional reduced basis spaces Xjv, Yv 
are constructed such that for any given \i £ T>, a solution (u*fi (fj,) , p 6 ^ (n)) £ Xn x Yv, 



f k {vN\ p), V V N £ X N , 



k £ 



(2.1) 



N ; 
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k e K, to (|2.1[) exists (see iJ2.ip , For /i € 2?, we then consider the errors 

e^(/i) = (e^ k (p,)) keK , where e^ k {p) = u £ ' k {^i) - u% k (fi) £ X, k € K, 

e p N (p) = (e^*(/i)) fc6K , where = p e ' fc (/z) - e Y, k e K, (2.3) 

e^(/i) = (e^ fe (^)) fceK , where e% (fi) = (e% (ji),e% (l*)) € Z, k € K, 

in the reduced basis approximations (u £ ^ k (fij^p 1 ^ (n)), k G K, with respect to the 
truth solution (u £,k (^,),p e,k (fi)), k g K; we note that in particular e^°(/i) = u e,0 (/i) — 
— from our initial conditions. 
To formulate our reduced basis a posteriori error bounds, we rely on the residuals 
associated with the reduced basis approximation (u e ^ (//) , p e ^ k (/*)) , fc 6 IK, 

- a(u% k (fi),v;n) - b(v,p £ jf(fi);n) € X', (2.4) 

r iv (sm) = 9 k ip) - K u n (m)><?;m) + £c(p S N k (v),q;v) e V, (2.5) 

for fc G K and n E T>. 

In the following analysis, we distinguish between the cases e = and e > 0. 

2.2.1. e = 0. We here derive rigorous upper bounds for the error e^(/x) measured 
in the "spatio-temporal" energy norm 

||(^) jeK ||/ a (o, fc; x)= ^ fe ||"+A«Ell uj H^J ' (4eiCI, fceK. (2.6) 

Our reduced basis a posteriori error bounds shall be formulated in terms of the dual 
norms of the residuals (|2.4p and (|2.5p . and (Online-)efScient lower and upper bounds 
to the truth continuity, coercivity, and inf-sup constants (jl.lOp . (II. lip . (I1.14p . and 



T^G") < 7m(*0 < 7^ B (M), 

7^(m) < 7a(M) < 7a UB (M), 

Q o B (a*) < OaOi) < a™(/i), 

/3 LB (^) < /3(M) < /? UB (m), 



V/ieP, (2.7) 



We can now state the following result. 

PROPOSITION 2.1. For any given fi £ V, N € N max , fceK, and a^ B (fi), 7™(a«) ; 
^ LB (m); 7m B (A*) satisfying (2l\) , we define 



A ^ lk^(-;M)II^ 2 / 7a B (^) \ I, i, jf ... „ aj, v„ 
At Z. ^b?^ +^b7^ U + ^bT^ ll^(^)IHI^(-;M)lk 



'^M , (7 a UB M) 2 Vi^(^)ii^ n l L 

At + a a B (m) J (/? LB (M)) 2 



(2.8) 



Then, Ajy(/x) represents an upper bound for the error e%(fi) measured in the "spatio- 
temporal" energy norm 12. 6\) . 

l|eft-M||*(o,fc;jc) < A^(/i), VfceK, ^eP, TVe N raax . (2.9) 
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Proof. Let \x be any parameter in T>, N <E N max , and k € K. For clarity of 
exposition, we suppress the argument /i in this proof. 

Take any 1 < j < k. From (23), ([23]) . and ([TTTSjl . the errors € X and 



AT 

i r- V cnfle-A, +v,„ + 

iV 



e^- 7 6 F satisfy the equations 



±m(e u J -e^-\v) + a(e^,v) + b(v,e p J) = r^(v), Vn£l, (2.10) 
By the LBB inf-sup condition ()1.17|) and (|2.10|) . we have 



b(e% j ,q)=r 2 J(q), V q e Y. (2.11) 



«ll pjii / Hv,eN ) r N (g) ~ Q ( e iV ; ") ~ Ai m ( e N ~ e N ; ") 
y < sup — — = sup — 

vex \\v\\x vex \\v\\x 
< \\r^\\x' + la\\e^\\x + ^\K j - etf~%, (2.12) 

where the last inequality follows from the Cauchy-Schwarz inequality for the inner 
product m(-, •), (TOT]) , and (1T"I"U|) . We then set v = e£ j , q = e p ^ in (pTTUl) . (j2"TTT) and 
subtract the second from the first equation such that 

J-m(p u ' j - p u ^~ X P u ' j ) 4- \\p u ' j \\ 2 — r 1 ^ (p u >°\ - r 2 J( p P>3\ 
At "H e ./V C N ' N I * \\ C N \\X,ii ~ 1 N \ C N I ' N I 



<\\r^\\x'\\e^\\ x + \\r^\\ Y ,\\e^\\ Y . 



Applying now (|2.12|) and (|1.14|l yields 

— m{e u ' J - - UJ_1 - UJ > —i— ii c u Jii 2 



^- ^ II IJll II 2,iii | (\\ . Ton 2,-7 II \ lle^llx,// 



1 V7m ||„2,j|| / 

At 



+ Hl r iV IMI e iV ~ C JV MMi 



which can be further bounded from Young's inequality by 



< i|l*lk-ll*llv. + ^(Wlk.+ * Hfbr-)' + ilk 



:Ll?IL|| r 2 ^|| 2 + _ i|2 

. 2 II AT 111" + 9A , 1 1 e 7V e W \\n- 



2At p 2 11 iV " r 2At 
Rearranging terms, the inequality now reads 



a. yi e JV lip 1 1 JV lljuj + W L N \\X,ii 

< !r» + if 1+ v\ i| r M|| x ,» r «i| v , + + ^ tMJt 



and the result follows from applying the sum Ylj=i> e N° = 0' an d (|2.7I) . □ 



8 



A.-L. Gerner and A. Reusken and K. Veroy 



In the special case of a symmetric problem, the error bounds given in Proposi- 
tion can be improved (see also [5]). We may then derive the following result. 

Proposition 2.2. Let a(-, ■; /i) be symmetric for all fj, G X>. For any given \x G T>, 
N G N max; k G K, and ^ B (p), 7? B (m)> /? LB (a4 7™ B (m) satisfying JO), we define 



a sym.fc / \ 



At 



k 

E 



"a B (f) 



/? LB (M) 




7 a UB (M) „ i. 



«a B (M) 



tT(m) 

At 



2-j 
'V 



(•;m)I 



(/3 lb (m)) 2 



1/2 



(2.13) 



Then, A s ^ m ' k (p) represents an upper bound for the error e%^(fi) measured in the 
"spatio-temporal" energy norm H2.6\) , 

(2.14) 



\\e u N ^)\\moMX) < tfT'Hn), V k G K, n G 2?, TV G 



Proof. Following the lines of the previous proof, we may now apply the Cauchy- 
Schwarz inequality for the inner product a(-, •) to obtain 

P\K j \\ Y < Wr^Wx- + v^ll^lk, + - 
instead of (|2.12p . Proceeding as before, this yields 

At m \ e N e N > e N ) + \\ e N lljf.M 

s II bill II 2,1'n , ( \\ r N ll-X' ' i V^ll 2,in \ II u,i|| 



1 



™ II r- 2 J II.,. II e " J _ p"'?" 1 ! 



At /3 

and the statement again follows from applying Young's inequality, the sum ; 
e u jf = 0, and $ZJ). D 

2.2.2. £ > 0. We here derive rigorous upper bounds for the error e e N (p) measured 
in the "spatio-temporal" energy norm 

/ k vVa 

||(« , ' s « , ')j g K|l^(o,fc i ^)= ^ fe ||"+A*Ell^ll^+ £ II^II^J . ( 2 - 15 ) 

where <7 J )jeK £ Z, fc G K. 

In addition to the dual norms of the residuals (|2.4p and (|2.5p , we here also rely on 
(Online-) efficient lower (and upper) bounds to the truth coercivity constants (jl.141) 
and (TTltl) . 



"a B W < «a(M) < «a B (M) 
C^) < ac(jU ) < a UB (M) 



lb,..w „ Jm : V ^e^ (2-16) 
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to formulate our reduced basis a posteriori error bounds. 

To demonstrate the differences to the case where e = 0, we may now recall the 
following result together with its proof (see [5]). 

PROPOSITION 2.3. For any given fi eV, N G N max , k G K, and a^ B (fi), a^ B (fi) 
satisfying \2.16\) . we define 

Then, A^ fe (/i) represents an upper bound for the error e £ N (/i) measured in the "spatio- 
temporal" energy norm &2.15\) , 

\\e e N (»)\\p(o,k;Z) ^ A '/(m), V k G K, [i G 2>, iV G N max . (2.18) 



Proof. Let /j, be any parameter in V, N G N max , and fc G K. For clarity of 
exposition, we suppress the argument /i in this proof. 

Take any 1 < j < k. From (|2~H) . and ([TTTgjl . the errors G X and 

G y satisfy the equations 

2n m ( e N 3 ~ e N 3 ~ 1 ' v ) + a ( e N J ^ v ) + ) = rl N( v ), V v £ X, 

H e N j iQ)- £C ( eP N>Q) = r N i (<l)> VqeF. 



Setting here v = e^ J , q = and subtracting the second from the first equation, we 
obtain 

i- t m(e# - e^" 1 , e#) + ||e^||^ + e ||e^||^ = r#(e#) - 

<l|r^||x'||^||x + K J '|k'|l4 J 'lk- ( 2 - 19 ) 
On the right-hand side, we now use ()1.14|) . (|1.15|) . and Young's inequality so that 

II^IU'||ey||x + ||r^||rHI^'l|y 

<MO^|| e ^|| +MpZi|| e w|| 

\ OL„. JOtr. 



< \ ( + ll^ll^ + + * lle^llU ; 

2 \ a a ea c J 

on the left-hand side, we use the Cauchy-Schwarz inequality for the inner product 
m(-, •) followed by Young's inequality so that 

™& d ~ e^~\e^) > \\e^\\l - \K j ~X\K% > \ (\KX " ll^ll*) ■ 
Rearranging terms, the inequality (|2 . 19[) finally reads 

and the statement follows from applying the sum Ylj=v = 0, and (|2.16l) . □ 
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2.3. Offline Online Computational Procedure. The efficiency of the re- 
duced basis method relies on an Offline-Online computational decomposition strat- 
egy. As it is by now standard, we shall only provide a brief summary at this point and 
refer to, e.g., [TTJ [24] for further details. The procedure requires that all involved op- 
erators can be affinely expanded with respect to the parameter [i. All /i-independent 
quantities are formed and stored within a computationally expensive Offline stage, 
which is performed only once and the cost of which depends on the large finite ele- 
ment dimension Af. For any given parameter [i £ V, the reduced basis approximation 
(u N k ([i),p N k ([j,)),k £ K, is then computed within a highly efficient Online stage; the 
cost does not depend on Af but only on the much smaller dimension of the reduced 
basis approximation space. The computation of the a posteriori error bounds con- 
sists of two components: The calculation of the residual dual norms \\r N (•; /i)||jf'> 
ll r jv fe ('i m)I|v > k G an d the calculation of the required lower and upper bounds (12.71) 
and (|2.16[) . respectively, to the involved constants. The former is again an application 
of now standard reduced basis techniques that can be found in [ill 121] • The latter is 
achieved by a successive constraint method (SCM) as proposed in [16]; we also refer 
to [7] for details in our saddle point context. 

2.4. Construction of Reduced Basis Approximation Spaces. We now 

turn to the construction of the reduced basis approximation spaces Xn,Yn, N £ 
N max . The low-dimensional spaces Xn,Yn are constructed by exploiting the para- 
metric structure of the problem: According to the so-called Lagrange approach, basis 
functions are essentially given by truth solutions associated with several chosen pa- 
rameter snapshots. However, in our time-dependent setting, Xn and Y/v not only 
have to appropriately represent the submanifold induced by the parametric depen- 
dence but also need to capture the causality associated with evolution in time to 
provide accurate approximations {u N k {p),p N k ([i)) for (u £ ' k ([i),p £ ' k ([j,)), k £ K, for 
any parameter query. Keeping computational cost to a minimum, we aim to achieve 
this with as few basis functions as possible. 

The POD greedy procedure represents an adaptive sampling process for parabolic 
problems that properly accounts for temporal and parametric causality: It combines 
the proper orthogonal decomposition (POD) method in k (see [2T] [22]) with the 
greedy procedure in \i (see [T] [3] and [7]). To begin with, we briefly recall the 
optimality property of the POD as described in [2T] [22]. For a given finite set 
Xx = {xi> ■ ■ ■ >Xz} Q A and Mx < dim(span(Ax)), the POD basis of rank Mx 
consists of Mx (■, -)x-orthonormal basis functions that approximate Xx best in the 
sense that 



analogously, we denote by PODy(3^r, My) the POD basis of rank My for a finite 
set y>x C Y, My < dim(span(3^r)). Assuming that we are given a current pair 
(Ajy_i, Y/v-i) of reduced basis approximation spaces, the POD greedy algorithm now 
proceeds as follows: In compliance with the greedy approach, it detects the parameter 
[in for which the (Online-)efficient reduced basis error bound attains its maximum 
over an exhaustive sample £ C T>. For a prescribed AJV £ K, we then compute the 
POD bases of rank AJV associated with the truth solutions u £,k ([iN) and p £ ' k ([iN), 
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k G K; more specifically, we compute PODjf (E u , AN) and POD Y (E p , AN) for 
E u = {u^ k ( f i N ) - U XN _ lU ^ k ( m ) | k G K}, 

^^{^M-ny-^p 6 '*^) I fcGK}, 

where and IIyv_ 1 refer to the (•, and (•, -)y-orthogonal projections on the 

current reduced basis approximation spaces Xjy-i and Y/v-i, respectively. Finally, 
the AX POD basis functions are appended to Xat-i and Y/v-i, respectively, and we 
obtain a subsequent pair (Xn,Yn). As before, this process is then repeated until a 
prescribed error tolerance is satisfied. We refer to [ITJ [HI [15] for a detailed dicussion 
of the POD greedy procedure, and to [HI [19] for an application to the Boussincsq 
and Fokker-Planck equations. 

For our saddle point problems, we now couple the above procedure with stabi- 
lization techniques developed in [7]; here (see also [9]), best convergence results were 
achieved by Algorithm 3 that aims to stabilize Xat, Yjy through an enrichment of the 
primal reduced basis approximation space with additional truth solutions. According 



Algorithm 1 Adaptive Sampling Procedure for e = 

l: Choose E c £>, dtoh *f c i G (0, 1), AN G K, and /zi G E 

2: Set TV 0, Vn <- {}, £>' <- {}, N Y <- 0, Fat <- {}, N x ^0,X N ^- {} 

3: repeat 

4: TV <- N+ 1, 2?jv "5- PjV-l U {/Zjv} 

5: EP = { P k { f X N )-U YN _ i p k ( t i N )\k€K} 

6: Ay <- Ay + AA, Fjv <- Fjv_i © span(PODy AA)) 

7: if fix then 

9: X x <- + AA, Aat <- JVat_i © span(POD x (£; tl , AA)) 

10: end if 

11: while (true) do 

12: for all fi G E do 

13: Compute (Ujy(jU),Pjv(jU)), fc G K, Ajv(a*), and 

14: (J*) = max { ^1-^^ ,0 } (cf. (23D)) 

15: end for 

16: ^ = argmax M6S A N (fi), {i* = argmax MeS d^(fi) 

17: if d$f(ji*) < <C then 

18: = Mat 

19: break 
20: end if 

21: if mn^evuvs > 0.1%, then 

22: V <- V U {/Z^} 

23: ^« = {7i fe ( A1 ' Ar )-n XivU fe ( A1 ' Ar )|fceiK} 

24: JV X f-JVx + AX, X N <- X N © span(POD x (£ u , AA)) 

25: else 

26: A A - N x + 1, X n <~ X n ® Bpaa{T ll .Q t r(jf) } (see (2.28), (2.36) in [7]) 

27: end if 

28: end while 

29: until A N (fj, N+1 ) < ^tol 

30: X max <- A 
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to these observations, we now apply the sampling procedures presented in Algorithm[T] 
and Algorithm [3J In case of e = 0, we use the distance d N (n) (see [7]), 

^Emj ^y }, ,ev, (2.20) 

of the inf-sup constants (3n(h) to the truth inf-sup constants /3(/i) as an indicator 
whether a current pair of reduced basis approximation spaces needs to be stabilized; 
the exact procedure is given in Algorithm [1] In case of e > 0, numerical results in 
[5] showed that the inf-sup constants /3/v(/z) may not be appropriate indicators for an 
ill-conditioned system but an adaptive sampling process should rather be based on 
the condition number k%(ji), 

e,max/ \ 

= ^ ■ W , e > 0, n G V, A G N max ; (2.21) 

here, cr^ max (/Lt) and Ojy min (/i) denote the maximum and minimum singular values 
of the corresponding reduced basis system matrix, respectively. Algorithm [2] now 
presents a possibility how this could be realized. 



Algorithm 2 Adaptive Sampling Procedure for e > 



Choose £ C V, 5 to \ G (0, 1), tf t re ol > 0, AN G K, and fi r G S 

Set TV <- 0, 2?x <- {}, 2?' <- {}, JVy <- 0, Y N <- {}, A x <- 0, ft <- {} 

repeat 

AT <- A + 1, V N <- Pat-i U {^jv} 

£f = { P £ ' fe (/ijv) - n^_ lP ^ fe (^) I k g k} 

Ay <- Ay + AN, Fat 4- Yx-i © span(PODy (JSP, AN)) 

if /ujv ^ 2?', then 

E» = {u e > k ( f i N )-IL Xtr _ 1 u e > k {n N ) I fcGK} 

Ax <— Nx + AiV, ft <- ATx-i © span(PODx(£", AJV)) 

end if 

while (true) do 
for all /i G £ do 

Compute (u^- (/i),^ (/*)), fc G K, Ajy(/x), and ^(/x) (see (|2.2ip ) 
end for 

Mjv = argmax MeS A n (/j,), [i* = argmax^sK^//) 
if K e N (fi* ) < <5« ol , then 

Mw+i = 
break 
end if 

if mm^v>uv N > 0.1%, then 

v <- v u 

= {^(^)-n Xjv u E ^(^) | fc g k} 

A x <- A x + A AT, X N ^ X N ® span(PODx(£ u , AN)) 
else 

Ax <^Nx + l,X N ^X N ® span{T^g N (fi*) } (see (2.28), (2.36) in [7]) 
end if 
end while 
until A N (fi N+1 ) < $tol 
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3. Model Problem. We consider a Stokes flow in a two-dimensional microchan- 
nel with an obstacle as introduced in [B] ; evolution in time is now induced by a time- 
dependent velocity profile on the inflow boundary. 

Let \x be any parameter in T>. For the physical domain Q and a given time interval 
[0, T] , T > 0, we now seek to find the (inhomogeneous) velocity M e ,inh : fi x (0, T) — > R 2 
and the pressure p c : £1 x (0, T) — > R satisfying 

- Au cMh + Vp c = in n x (0, T), (3.1) 

V-u e ,inh = infix (0,T), (3.2) 
subject to initial conditions u Cj i n h(-,0) = and with boundary conditions 
u e ,mh(£, t) = H(t)h(x) on T in x (0, T), u e , inh = on f o x (0, T), 



9h 



p e n on r out x (0,T); 



here, A and V again denote the Laplacian and gradient operator over the physical 
domain (l, n is the unit outward normal, h : R 2 — > R 2 is given by h(x) = (4x2(1 — 
x 2 ), 0) for all x = ( Xl ,x 2 ) £ R 2 , and we choose H : [0, T] -> R with ff(t) = t(sin(27rf)+ 
1) for all t £ [0,T]. According to the setting introduced in [6], we also consider the 
following perturbation of the problem (|3.1j) - (|3.3j) : For a sufficiently small e > 0, we 
again introduce a penalty term to the continuity equation (|3.2p such that 

V-ul Anh = -s P l inQx(0,T). (3.4) 

We now follow the steps discussed in [B]: We choose the lifting function u? = 
Hul where ul is defined as in [6J, and transform the problem statement for the 
homogeneous velocity u E c = u e c inh — to an equivalent problem posed over the 
reference domain tt. Furthermore, as required for the time-discrete setting introduced 
in fJU we divide the time interval [0, T] into K subintervals of equal length At = T/K, 
and consider a backward Euler method for time integration. The problems (|3.1I) - (|3.3I) 
and (|3.ip . (|3.4p . (|3 . 3[) may thus be written as a parametrized saddle point problem of 
the form (|1.9[) . Here, for any /i £T>, the bilinear forms a(-, •; /i), &(•, •; /i), and c(-, •; /i) 
are given as in [BJ; accordingly, the bilinear form m(-, •; /i) : X c x X — > R represents 
the L 2 -inner product for vector functions over the physical domain Q, formulated on 
the reference domain J7, 

and the linear functionals /(•; fi) and g(-; fj.) are given by 

f(v,t; t i) = f(v,t) = -H'(t)[ u h -vdx-H(t)[ ^dp-dx, 

Jn L oxj ox j 

9 (<1, *; A*) = 9 (?)*)= H(t) [ q^^-dx, 

Jo L oxi 

for all v £ X c , q £ Y c , t £ [0,T]. We recall that the bilinear forms a(-, ■; /i), &(■, •; /i), 
and c(-, •; /i) then satisfy the assumptions (|1.2I) - (|1.8I) . (|1.4I) . and (|1.6|) . For all /ieD. 
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m(-. defines an inner product on X e such that (|1.7p holds true; moreover, there 
exists a constant C|(/i) > from the Poincare inequality (see, e.g., (53]) such that 

m(v, v; fi) < C%{p) a(v, v; fi), V v e X c , V fi G V, 

and thus (jl.ip is satisfied with 

7m (/•*) = sup sup ' ' — < Q(m)72G") < °°> V neV. 
uex c vex \\u\\x e \\v\\ Xe 

We particularly have Q a = 10, Qb = 6, Qf = Qg = 1, and Q m = Q c = 5 in the 
respective /z-affine expansions (see £12.31) . 

Choosing the truth approximation spaces X and Y as the standard conforming 
P2-P1 Taylor-Hood finite element approximation subspaces [25] over the regular tri- 
angulation 7a, we ensure that also (11.171) is satisfied (see, e.g., [2] [4] [lOj [23] ) and 
therefore recover the situation described in §1.21 

4. Numerical Results. We now apply the reduced basis methodology devel- 
oped in <J5] to our model problem introduced in $3] We set T — 1 and consider a 
constant time step size At corresponding to K = 100 time levels. The truth discretiza- 
tion is based on a fine mesh with a total of J\f = 72,076 velocity and pressure degrees 
of freedom. In this section, all numerical results are attained using the open source 
software rbOOmit [20j . an implementation of the reduced basis framework within the 
C++ parallel finite element library libMesh [17] . 

4.1. e — 0. We first turn to the coercivity, continuity, and inf-sup constants 
required for our reduced basis procedure. We obtain (Online-)efficient lower and upper 
bounds to a a (fx), -fain), and by using the SCM (see §2.3|) with the configurations 
specified in [TJ . To estimate the continuity constants j m (//) , we apply the method for 
M a — 00, M + — 0, an exhaustive sample 5 C T> of size |S| = 4,225, and the SCM 
tolerance e = 0.01 (see [H]). We then obtain accurate (Online-) efficient lower and 
upper bounds 7 m B (M) and 7 m B (£i) with if max = 5. 

We now turn to the reduced basis approximation. To build our low-dimensional 
reduced basis approximation spaces Xn,Yn, N 6 N max , we apply the POD greedy 
procedure described in Algorithm[T](see £|2.4[) . The sampling process is based on an ex- 
haustive random sample E C T> of size |E| = 4,900, AiV = 2, and S^ ol — 0.1; since our 
Stokes model problem is clearly symmetric, wc here in particular use the relative re- 
duced basis a posteriori error bound Ajv(a*) = A s ^ m ' K (u' N (fi))jg&\\£2(o t }c ] x) (see 
(12. 6p . (|2.13p ). Figure |4~T1 now shows the maximum error || e]^- (/i) ||^2(o,ftT; jc) (see (I2.3p ) 
in the reduced basis velocity approximations and associated error bounds A s ^ m ' K (fx) 
and A^(/i) (see (|2.8p ) as functions of the dimension Nz\ Figure l4~2l presents the 
maximum error ||e^ r (/i)||^2( ,fe;X) and associated error bounds A s ^ m ' k (fi), A^,(/i) as 
functions of k £ K for several values of Nz ■ 

First, we observe that the reduced basis error and error bounds are roughly uni- 
form in time (see Fig. 14. 2p and decrease rapidly as Nz increases (see Fig. 14. ip . We 
obtain stable, rapidly convergent reduced basis approximations, and rigorous a poste- 
riori error bounds that reflect the behavior of the error very accurately. Second, the 
error bounds are tight. To quantify this statement, we present in Table 14+1 maximum 
effectivities associated with A^ m ' fc (/i) and A^(/i) for several values of k and N. We 
notice that their values remain more or less constant with k. Moreover, as in the sta- 
tionary case (see [5]), we again benefit from exploiting the symmetry of the problem: 
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Fig. 4.1. Maximum error He^y-dOIL^ro kx) ( see {E2P j {SUB) an d maximum error bounds 

A^ m ' (/i) and A^(^t) ( see 12.1 tit and 12. 8j ) normalized with respect to |]('U-'(/u))jeKllf 2 (o,X;X) 
shown as functions of Nz; the maximum is taken over 25 parameter values. 



N = 13 {N z = 76) 



l|e"ll 
. . . \ s v m -' 



10 20 30 40 50 60 70 

k 



AT = 20 (JV Z = 109) 



N = 41 (JV Z = 226) 



10 20 30 40 50 60 70 80 90 100 10 20 30 40 50 60 70 80 90 100 

k k 



Fig. 4.2. Maximum error ||e^ r (A t )||^2(g k-X) ( see i2.3t , 12. 6V ) and maximum error bounds 
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(/i) anrf A^(/j) (see \2. 13V and i2.8V ) normalized with respect to || (it- 5 ' (A*))jeKl 



l 2 (0,k;X) 



shown 



as functions of k € K for several values of N; the maximum is taken over 25 parameter values. 



Effectivities range from 33 to 45 in case of A^O) (see TableO (b)) and improve m 
case of A^ m ' fe (yLt) by roughly 10 (see Table |4~T1 (a)). We emphasize at this point that 
the error bound formulations in (I2.8[) and (I2.13[) in fact suggest a growth in time. In 
practice, this behavior could not yet been observed and may therefore be investigated 
in greater detail within future work. 

We now discuss the Online computation times for the proposed method. For 
comparison, once the /x-independent parts in the affine expansions of the involved op- 
erators have been formed (see $2.3\i , direct computation of the truth approximation 
{u k (/i) , p k (/i)) , k € K, (i.e., assembly and solution of (|1.18[) ) requires roughly 30 sec- 
onds on a 2.66 GHz Intel Core 2 Duo processor. We initially take a total reduced basis 
dimension of Nz = 226. Once the database has been loaded, the Online calculation of 
(itjy(/z),pjy(/z)), k € K, (i.e., assembly and solution of (|2.1|l ) and A s ^ m ' k (fi),k £ K, for 
any new value of \i £ T> takes on average 27.97 and 80.76 milliseconds, respectively, 
which is in total roughly 270x faster than direct computation of the truth approxima- 
tion. Thus, even for this large value of Nz, we obtain significant Online savings. In 
practice, however, wc quite often need not take such a large value of Nz', our rigorous 
and inexpensive error bounds A^ m ' fc (/x), k £ K, allow us to choose the reduced basis 
dimension just large enough to obtain a desired accuracy. To achieve a prescribed ac- 
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(a) Effectivities rfS m ' (fJ,) associated with A^X m ' (/j) 



N 


N z 


k = 10 


fc = 20 


fc = 40 


fc = 60 


fc = 80 


fc = 100 


5 


36 


28.90 


30.38 


31.30 


31.61 


31.63 


31.13 


10 


63 


30.41 


31.35 


32.20 


32.37 


32.16 


31.91 


15 


93 


25.83 


28.05 


29.39 


29.51 


29.44 


29.38 


20 


109 


23.31 


24.25 


26.79 


27.23 


27.08 


27.21 


25 


142 


25.29 


28.15 


29.54 


29.73 


29.66 


29.64 


30 


177 


26.28 


26.05 


28.77 


30.58 


30.70 


30.60 


35 


201 


24.77 


24.86 


26.68 


27.36 


27.51 


27.81 


40 


222 


24.18 


23.96 


24.19 


25.03 


25.51 


25.54 


(b) Effectivities r]% (/J.) associated with A^(/t) 
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N z 
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5 


36 
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42.76 


43.20 
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42.53 


10 


63 
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41.66 
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43.00 


42.90 


15 


93 


38.58 
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45.14 


45.07 


20 


109 


32.59 


34.10 


37.11 


37.60 


37.17 


37.19 


25 


142 


35.52 


39.26 


42.93 


42.73 


42.63 


43.55 


30 


177 


34.31 


34.27 


36.41 


37.31 


39.67 


40.64 


35 


201 


32.86 


33.59 


36.52 


37.39 


36.65 


37.54 


40 


222 


33.67 


33.76 


35.25 


35.42 


34.78 


35.17 



Table 4.1 

Maximum effectivities (a) n B ^ In ' k (fi) = A^" 1 '* (fJ.)/\\e^ (/x) \\ e 2 (q u-X) ( see t& 14t J an d 

(b) v^(fi) = A^(/^)/||e^r(/x)|| f 2(o k . x j (see I2.9\l ) for several values of fc S K and N; the maxi- 
mum is taken over 25 parameter values. 



curacy of at least 1% (resp., 0.1%) in the reduced basis approximations u k N {^), k G K, 
we need Nz = 76 (resp., Nz = 109) (see Fig. 14. ip . Again, once the database has 
been loaded, the Online calculation of fc G K, and A^ m ' fe (fi) , k G K, 

for any new value of \i G V takes on average 4.41 (resp., 7.62) and 24.47 (resp., 33.75) 
milliseconds, respectively, which is in total roughly 1,000 x (resp., 700 x) faster than 
direct computation of the truth approximation. 

4.2. e > 0. Again, the SCM (see §2.3|) enables the (Online-)efficient estimation 
of the coercivity constants a a {fi) and a c (/n); as we here use the same configurations, 
we refer to [5] for details in this context. 

To build our low-dimensional reduced basis approximation spaces Xjf, Yjv, TV e 
N max , we apply the POD greedy procedure described in Algorithm [2] (see $2A\i . The 
sampling process is based on an exhaustive random sample S C T> of size |S| = 
4,900, A7V = 2, 6£ ol — 10 3 , and the relative reduced basis a posteriori error bound 

AjvOu) = ^% K (p)/\\(u e ^{/i))jOL\\p(o,K-,Z) (see ^M, (233). Figure Efl now shows 
the maximum error ||e|/(A*)||^(o ,K;Z) ( se e (|2.3[) ) in the reduced basis velocity and pres- 
sure approximations together with the associated error bound (/i) as functions 
of the dimension Nz for different values of e. Figure l4~4l then presents the maximum 
error ||e|^(/i)||^2( 0> fc ; z) and associated error bound A^ fc (/z) as functions of fc G K for 
several values of N; note that the latter arc chosen as the values for which the er- 
ror bounds A e ^ K (fi) guarantee a prescribed accuracy of at least 1% and 0.1% in the 
reduced basis approximations. 

First, we again observe that the reduced basis error and error bounds are roughly 
uniform in time (see Fig. I4.4p and decrease rapidly as Nz increases (see Fig. I4.3|) . We 
obtain stable reduced basis approximations whose rapid convergence is not affected by 
the penalty parameter, and a posteriori error bounds that are meaningful and rigorous. 
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Second, using the condition numbers K s N (fJ-) as an indicator for an ill-conditioned 
system, Algorithm [2] guarantees stability by properly accounting for the effects of 
the penalty term: For e = 10~ 2 , the sampling process recognizes that the reduced 
basis approximation spaces Xn , X/v do not have to be stabilized to provide accurate 
approximations; taking smaller values of e and thus approaching the non-penalized 
problem, an additional enrichment of the reduced basis approximation space for the 
velocity then becomes more and more necessary. Third, we see that the error bounds 
are tight for e = 10 -2 but become less sharp as we decrease e and our perturbed 
truth approximation becomes more accurate. However, effectivities exhibit a similar 
O(-^j) -dependence on the penalty parameter as observed in the stationary case (see 
[B]) and remain reasonably small for relatively small values of e. To further quantify 
this statement, we present in Table POl the effectivities associated with A^ fc (/i) for 
different values of k, N, and e. We note that their values are fairly constant with 
k and N and confirm the O(-^j) -dependence indicated by Fig. 14.31 and Fig. 14.41 As 
before, the effects of the penalty parameter on the effectivities are thus relatively 
benign and we obtain useful bounds for reasonably small values of e. 

We close this section by discussing the Online computation times. For com- 



e = icr 2 e = l(r 3 




N z 



Fig. 4.3. Maximum error \[e.%[pi)\\i^iQ t K;Z) ( see lEM i \SHM) an d maximum error bound 
A £ ^ K (fi) (see \2.17\ l) normalized with respect to \\ (u e ' J (fi), p e ' J '(m))j'6K \\e 2 (o K-Z) shown as functions 
of Nz for different values of e; the maximum is taken over 25 parameter values. 
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parison, once the //-independent parts in the affine expansions of the involved op- 
erators have been formed (see b|2.3|) . direct computation of the truth approximation 
(u £ ' k (fi),p £ ' k (/i)),k £ K, (i.e., assembly and solution of (|1.18[1 ) requires roughly 23 
seconds on a 2.66 GHz Intel Core 2 Duo processor. Again, our rigorous and inex- 
pensive reduced basis a posteriori error bounds enable us to choose the reduced basis 
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Fig. 4.4. Maximum error |[eJ^(/Lt)[|^2/Q t k;Z) ( see E2P> f2.150 ) and maximum error bound 
A £ ^- k (fj,) (see $2.17]) ) normalized with respect to \\ (u E ' J (fJ.), p e ' J ; (M))j'eK \\i 2(o k-z) shown as functions 
of k G K for several values of N for (a) e = 10~ 2 , (b) e = 10" 3 , (c) e = 10~ 4 , and (d) e = 10~ 5 ; 
the maximum is taken over 25 parameter values. 
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dimension just large enough to obtain a desired accuracy. Choosing e = 10~ 2 , the er- 
ror bounds A £ j^ k (fj,) are sharp with effectivities of approximately 12 (see Table [PI (a)) 
and prescribe a dimension of Nz = 68 to achieve an accuracy of at least 1% in the re- 
duced basis approximations (u £ ^ k (/j,) , p £ ^ k (fi)) , k e K (see Fig. l4.3|) . Once the database 
has been loaded, the Online calculation of (u 5 ^ ([i) , pjf (/j,)) , k £ IK, (i.e., assembly and 
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Table 4.2 

Maximum effectivities rjjj (jjl) = A^ fc (/i)/||e| sr (/i)||^2(Q fc; ^) (see j2. 18)) ) for several values of 
k £ K and N for (a) e = 10 — 2 , (b) s = 10 -3 , (c) e = 10 — 4 , and (d) e = 10 — 5 ; the maximum is 
taken over 25 parameter values. 
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solution of (|2.ip ) and A^r (/x), fc £ K, for any new value of € V then takes on average 
3.71 and 14.53 milliseconds, respectively, which is in total roughly l,200x faster than 
direct computation of the truth approximation. Choosing smaller values for e, the 
error bounds become more pessimistic and thus dictate a larger system dimension at 
which they guarantee the same order of accuracy. For e = 10 -5 , we need Nz = 121 to 
achieve a prescribed accuracy of at least 1% in the reduced basis approximations (see 
Fig. 133]); the Online calculation of (u^ (fi) , p e ^ h (//)) , fc e IK, and A^ fe (/i), fc e IK, then 
takes on average 9.19 and 33.81 milliseconds, respectively, which is in total roughly 
500 x faster than direct computation of the truth approximation. Thus, even for small 
penalty parameters e, accurate approximations are guaranteed at significant Online 
savings. Detailed computation times for different values of e are given in Table 14.31 
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Table 4.3 

Average computation times in milliseconds for the Online evaluation of (ji^ fe (/x), p^ fe (^t)), fc 6 
K, (assembly and solution of 12. U ) ) and the error bounds A^ k ([i), fc £ K, (see \2.1T\j ) for different 
values of e with a prescribed accuracy of at least 1% (resp., 0.1%) for the reduced basis approxima- 
tions («^ fc (^t),p^ fe (/i)), tel. 



5. Concluding Remarks. In this paper, we present new reduced basis methods 
for the instationary Stokes equations. 

Combining techniques developed in [7J |S] with current reduced basis approaches 
for parabolic problems, we derive new rigorous a posteriori bounds for the errors in 
the reduced basis velocity approximations and a POD greedy procedure that properly 
accounts for temporal and parametric causality as well as stability. The method 
provides rapidly convergent reduced basis approximations that are highly efficient 
and whose accuracy is certified by sharp and inexpensive a posteriori error bounds. 

An approximation by penalty or regularization again allows for significant Offline 
savings at the expense of a less accurate truth approximation. Due to the introduced 
penalty term, an additional enrichment of the reduced basis velocity approximation 
space is not always necessary to obtain stable approximations; moreover, we obtain 
a posteriori error bounds that do not involve the expensive computation of inf-sup 
stability constants. As in the stationary case (see [B]), the method provides reduced 
basis approximations and meaningful a posteriori error bounds that are computed 
very easily; nevertheless, drawbacks such as the disadvantageous dependence of the 
error bounds on the penalty parameter remain. 

Time integration is achieved through a backward Euler method. Clearly, also 
other time integration schemes may be used. Using a Crank-Nicolson method, often 
preferred in practice due to its second-order accuracy, we may develop a penalty 
approach that is very similar to the one presented in this paper (see [5]); in case of 
£ = 0, useful reduced basis a posteriori error bounds could not yet been derived and 
may therefore be — as well as a posteriori error bounds for the reduced basis pressure 
approximations — part of future work. 
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